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We search for physics beyond the standard model using events with a photon, two or more 
hadronic jets, and an apparent imbalance in transverse energy, in pp collisions at the Fermilab Teva- 
tron at \fs = 1.8 TeV. Such events are predicted for production of supersymmetric particles. No 
excess is observed beyond expected background. For the parameter space of the minimal supersym¬ 
metric standard model with branching fraction B(x§ —> 7Xi) = 1 and rn^o — m^o > 20 GeV, we 
obtain a 95% confidence level lower limit of 310 GeV for the masses of squarks and gluinos, where 
their masses are assumed equal. 


2 



We search for physics beyond the standard model (SM) 
using events with one high transverse energy ( Et ) pho¬ 
ton, two or more jets, and large imbalance in transverse 
energy We call these + > 2 jets events. This 

search is motivated by recent suggestions |l],|2j] that su¬ 
persymmetry may result in signatures involving one or 
more photons together with multiple jets and large Ipr- 

Supersymmetry is a generalization of space-time sym¬ 
metry. It introduces for every particle in the standard 
model a supersymmetric partner differing in spin by 
one half. A-parity Q, defined as +1 for SM particles 
and —1 for their super-partners, is assumed to be con¬ 
served in this analysis, such that supersymmetric parti¬ 
cles are produced in pairs and the lightest supersymmet¬ 
ric particle (LSP) is stable. In the minimal supersym¬ 
metric standard model (MSSM), the gaugino-Higgsino 
sector (excluding gluinos) is described by four param¬ 
eters: Mi, M 2 , p, and tan/3, where Mi and M 2 are 
the U(l) and SU(2) gaugino mass parameters, p is the 
Higgsino mass parameter, and tan/3 is the ratio of the 
vacuum expectation values of the two Higgs doublets. 
Gaugino-Higgsino mixing gives four neutral mass eigen¬ 
states (neutralinos x°, i = 1, ...,4) and two charged mass 
eigenstates (charginos xf , i = 1,2). Within the MSSM, 
the radiative decay of x 2 7Xi dominates in the re¬ 
gion 50 < Mi ~ M 2 5 100 GeV, 1 < tan/3 < 3 and 
—65 ;$ p —35 GeV of parameter space Q, and has 
been proposed as an explanation Q of a candidate event 
reported by the CDF Collaboration ||. Assuming that 
Xi is the LSP, then the production of x 2 , either directly 
or from decays of other supersymmetric particles, will 
yield 7 J$t + -V events. 

In this Letter, we present a search for physics beyond 
the SM in the channel pp — > 7 IpT + > 2 jets at the Fer- 
milab Tevatron collider. Because of large backgrounds 
from QCD processes, we do not consider events with less 
than two jets. We interpret our results in terms of squark 
(q) and gluino (g) production in the context of supersym¬ 
metric models with a dominant x 2 —> 7 X 1 decay. 

The data used in this analysis were collected with the 
D0 detector during the 1992-1996 Tevatron run at a 
center of mass energy of y/s = 1.8 TeV, and represent 
an integrated luminosity of 99.4 ± 5.4 pb _1 . A detailed 
description of the D0 detector can be found in Ref. ||. 
The trigger requires one electromagnetic (EM) cluster 
with Et > 15 GeV, one jet with Et > 10 GeV, and 
]$t > 14 GeV {Ifir > 10 GeV for about 10% of the data 
taken early in the Tevatron run). Photons are identified 
via a two-step process: the selection of isolated EM en¬ 
ergy clusters, and the rejection of such clusters with any 
associated charged tracks. The EM clusters are selected 
from calorimeter energy clusters by requiring: (i) at least 
95% of the energy to be deposited in the EM section 
of the calorimeter; (ii) the transverse and longitudinal 
shower profiles to be consistent with those expected for 
an EM shower; and (iii) the energy in an annular isola¬ 


tion cone with radius ( 1Z = 1 J (A(/>) 2 + (A?/) 2 ) 0.2 to 
0.4 around the cluster in p — cj> space to be less than 10% 
of the EM energy in an 1Z = 0.2 cone, where rj and (j) 
are the pseudorapidity and azimuth, respectively. The 
EM clusters that have either a reconstructed track or 
a large number of hits in the tracking chamber along a 
road joining the cluster and the interaction vertex are 
vetoed. I£t is determined from the energy deposition in 
the calorimeter within |^| < 4.5. 

To be selected as xtfr + > 2 jets candidates, events 
are first required to have at least one identified photon 
with Et > 20 GeV and pseudorapidity |?/ 7 | < 1.1 or 
1.5 < |r/ 7 | < 2.0, and two or more jets reconstructed with 
cones of radius 1Z = 0.5, having E 3 T > 20 GeV and rp < 
2.0. We refer to the events passing these requirements as 
the 7 + > 2 jets sample. The ]$t distribution of these 
events is shown in Fig. |l]. We then require Ifr > 25 GeV. 
A total of 318 events satisfy all requirements. 



FIG. 1. The Jpr distributions of the 7 + > 2 jets 
(solid circles) and background (solid histogram) events. The 
number of events in the background is normalized to the 
7 + > 2 jets sample for Ifr < 20 GeV, the region left of 
the dot-dashed line. Also shown (dashed and dotted his¬ 
tograms) are the distributions expected from supersymmetry 
for rriq = nig = 150 GeV and 300 GeV. 

The principal backgrounds to the signal are: events 
from sources such as QCD direct photon and multijet 
events, where there is mismeasured I^t and a real or fake 
photon; W(—> ev) + jets events, where the electron is 
misidentihed as a photon; and W (—> £v) + jets events 
(where £ = e, /r, r), in which one of the jets is misidenti- 
fied as a photon. These backgrounds are estimated from 
the data sample with the same trigger as the candidate 
events. The background from mismeasurement of J$t 
is estimated using events with one EM-like cluster that 
satisfies all photon criteria, except requirement (ii) on 
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the shower profile. These events must also have two or 
more jets with E 3 T > 20 GeV and \r/ 3 \ < 2.0, making 
them similar to those of the 7 + > 2 jets sample, and 
therefore of similar resolution in IpT- The events in this 
background sample are normalized to the 7 + > 2 jets 
sample for ]fir < 20 GeV, which provides an estimated 
background from mismeasurement of 315 ± 30 events 
beyond 1$t = 25 GeV. 

W + > 2 jets events with W —> ev can mimic 7 Ifir + > 
2 jets events if the electron is misidentified as a photon. 
This contribution is estimated using a sample of elpT + > 

2 jets events that passes all our kinematic requirements, 
with the electron satisfying those defined for the photon. 
Electrons are selected from identified EM clusters that 
have matched tracks. The probability that an electron 
is misidentified as a photon is determined from Z —> ee 
events as 0.0045±0.0008. Multiplying this probability by 
the number of etfr + > 2 jets events yields a background 
of 4 ± 1 events. 

The W (—> Iv) + jets background is estimated using a 
data sample of W {—> ev) + > 3 jets events passing all 
kinematic requirements, with at least one of the jets sat¬ 
isfying those imposed on photons. The probability that a 
jet is misidentified as a photon is determined by counting 
the number of photons observed in multijet events. We 
find this to be 0.0007±0.0002. Using this probability and 
the scale factor N W (_> tu ) + > 3 jets /N w ^ ev)+ > 3 jets (deter¬ 
mined from Monte Carlo), we estimate a background of 
1.0 ± 0.3 events. The background from Z{— > vv) + > 

3 jets is found to be negligible. 


Number 

IpT> 25 GeV 

25 GeV 

I$t> 50 GeV 

of 

No Ht cut 

Ht > 200 GeV 

No H t 

cut 

jets 

N S 

Nb 

N S 

Nb 

N S 

Nb 

n > 2 

318 

320±30 

30 

20±10 

43 

65±15 

n > 3 

70 

70±15 

17 

8±5 

11 

10±5 

n > 4 

8 

10±5 

6 

4±3 

1 

3±3 


TABLE I. Number of observed 7 Ifx + n jets events (Ns) 
together with the corresponding number of background events 
(Nb) for n > 2,3,4, for three sets of cutoffs. 


The number of observed events and the expected back¬ 
grounds are summarized in Table ||, together with break¬ 
downs into events with three or more and four or more 
jets. The Ht distribution (defined as the scalar sum of 
the Et of all jets with E° T > 20 GeV and rf < 2.0) is 
shown in Fig. [|, for both 7 Ifix + > 2 jets and background 
samples. The background distribution is consistent with 
that observed for 7 J$t + > 2 jets. Also given in Table | 
is the number of observed events and the expected back¬ 
ground if the cutoff Ht > 200 GeV is applied or if the 
I^t cutoff is raised to 50 GeV. In all three comparisons, 
the estimated number of background events agrees with 
the number of events observed in the data. 

To optimize selection criteria for a supersymmetric sig¬ 
nal, we simulate squark and gluino pair production, and 
also production in association with charginos or neutrali- 



FIG. 2. The Ht (defined as y . E 3 T ) distributions of the 
7 IpT + > 2 jets and background events. The expected distri¬ 
butions from supersymmetry are also shown for comparison. 

nos using the SPYTHIA program 0. The MSSM pa¬ 
rameters are set to M\ = M 2 = 60 GeV, tan/3 = 2, 
and n = —40 GeV. This set gives m^o = 34 GeV, 

rri^o = 60 GeV, and B(%2 —> 7X1) = 1- Sleptons (I) 
and stop (t\) are assumed to be heavy. Events with x 2 
in the final state are selected and processed through the 
D0 detector-simulation program and the trigger sim¬ 
ulator. The same trigger requirements, reconstruction, 
and selection criteria are then applied as were used with 
the data. Monte Carlo (MC) events are generated for 
three squark or gluino mass possibilities: (i) equal mass 
squark and gluino ( rriq = rrig); (ii) heavy squark and 
light gluino (rriq rrig ); and (iii) light squark and heavy 
gluino (rriq <C Trig). The Ifir and Ht distributions for 
m.q = m 3 = 150, 300 GeV events are shown, respectively, 
in Figs, jl] and ||, where the MC distributions are scaled 
by the factors shown in parentheses. The distributions 
expected from supersymmetry differ considerably from 
those of the background. To increase the sensitivity to 
supersymmetry, we introduce an Ht cutoff, and maxi¬ 
mize the 6s/SNb ratio by varying the IpT and Ht cutoffs. 
Here es is the efficiency for signal, and 6Nb is the un¬ 
certainty on the estimated number of background events. 
To ensure high efficiencies for both low and high squark 
and gluino masses, the optimization is done for two MC 
points rriq = rrig = 150 and 300 GeV. The optimum val¬ 
ues are 1^t> 35 GeV and Ht > 100 GeV for 150 GeV, 
and #r> 45 GeV and H T > 220 GeV for 300 GeV. 
The es/SNb results (a function of squark/gluino mass) 
for the two sets of optimized cutoffs are equal near 
200 GeV. Therefore we apply the cutoffs optimized for 
the 150 GeV mass point to MC events with squark and 
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gluino masses below 200 GeV, and apply those optimized 
for the 300 GeV mass point to masses of 200 GeV or 
above. The number of events observed for these two sets 
of cutoffs are 60 and 5, with 75 ± 17 and 8 ± 6 events 
expected from background processes. We consequently 
observe no excess beyond the standard model. 


rn g/9 

GeV 

mq(= 

e 0 (%) 

mg) 

e S (%) 

mg(<g mq) 

eo (%) e s (%) 

m 5 -(C mg) 

eo (%) e s (%) 

150 

66.2 

15.1i0.8 

69.1 

11.6T0.9 

60.0 

16.8il.l 

200 

62.3 

7.9i0.6 

59.6 

5.34=0.6 

53.8 

9.5i0.9 

250 

59.6 

14.8T0.8 

49.7 

13.6=1=1.1 

55.4 

14.8il.l 

300 

56.1 

21.5il.O 

43.1 

19.Oil.3 

55.4 

22.1T1.2 

350 

51.8 

22.8±1.1 

39.3 

23.5T1.5 

52.7 

26.6il.4 

400 

46.7 

23.5T1.1 

33.3 

22.7±1.6 

54.3 

25.8il.3 


TABLE II. The percentages of events (eo) generated con¬ 
taining x° in the final state, and the efficiencies (es) for their 


detection using the two sets of optimized cutoff's as discussed 
in the text, for different values of squark/gluino mass. The 
uncertainties are purely statistical. 


The detection efficiencies (eg) for predictions from the 
supersymmetric models are given in Table |J, along with 
the percentages (eo) of generated events having \2 i n 
the final state. MC studies show that the overall ef¬ 
ficiency varies by 4% for different choices of Mi, M 2 , 
tan/3, and p that are consistent with B(x 2 —> 7 X?) = 1 
and to^o — m^o > 20 GeV, the suggestions offered in Ref. 
Q. Experimentally, the mass requirement is needed to 
ensure that photons from \2 decays are reasonably en¬ 
ergetic and can be detected with good efficiency. The 
total systematic error on the efficiency is 9%, including 
uncertainties in photon identification efficiency (7%), the 
choice of values of the supersymmetry parameters (4%), 
and the jet energy scale (3%). 

We set a 95% confidence level (C.L.) upper limit on 
cr x B = cr(pp -> q/g -> x° + A') x B(x° -> 7 %?) us¬ 
ing a Bayesian approach with a flat prior distribution for 
the signal cross section. The statistical and systematic 
uncertainties on the efficiency, the integrated luminosity, 
and the background estimate are included in the calcula¬ 
tion of the limit, assuming Gaussian prior distributions. 
The resulting upper limit as a function of squark/gluino 


mass is tabulated in Table III 


m g/<7 

GeV 

mg(= 

Theory 

mg) 

Limit 

cr X B (pb) 

mg) 

Theory Limit 

mg(<Si 

Theory 

mg) 

Limit 

150 

83.4 

2.0 

24.1 

2.6 

8.51 

1.8 

200 

12.1 

1.1 

3.48 

1.6 

1.59 

0.9 

250 

2.37 

0.57 

0.51 

0.63 

0.43 

0.58 

300 

0.53 

0.39 

0.12 

0.44 

0.12 

0.38 

350 

0.13 

0.37 

0.02 

0.37 

0.03 

0.32 

400 

0.04 

0.36 

0.01 

0.37 

0.01 

0.32 


TABLE III. The theoretical cross section a x B and our 
measured 95% confidence level upper limit onaxB for differ¬ 
ent values of squark/gluino mass. The predictions are calcu¬ 
lated for Mi = M 2 = 60 GeV, tan /3 = 2, and p = —40 GeV. 


Figure shows the limit for the case where rriq = nig, 


together with the leading-order theoretical cross section, 
calculated using the SPYTHIA program with the CTEQ3L 
parton distribution functions |J. The renormalization 
scale (prs) is set to the average transverse energy (( Et)) 
of the outgoing partons in the calculation. The cross 
section varies by about ±30% if pus = 2 (Et) or pus = 
(Et)/2 is used. The hatched band represents the range 
of predictions obtained by varying the supersymmetry 
parameters with the constraints that B(x° —> 7 X?) = 1 
and m^o — m^o > 20 GeV, assuming prs = (Et)- The 
intersection of the limit with the lower edge of the band 
is at cr x B = 0.38 pb, leading to a lower limit for equal 
mass squarks and gluinos of 310 GeV at the 95% C.L. 



m(q) (GeV) 


FIG. 3. The 95% C.L. upper limit on a x B as a function of 
m,g/g, assuming equal squark and gluino masses. The hatched 
band represents the range of expected cross sections for differ¬ 
ent sets of MSSM parameters, consistent with the constraints 
B(x§ —> 7 X 1 ) = 1 an d rn^o — m^o > 20 GeV. The inflection 
below 200 GeV in the limit curve is the intersection of the 
two curves using the two sets of optimized cutoffs discussed 
in the text. 

The effect of light sleptons on squark and gluino de¬ 
cays is studied by varying the slepton mass (to| = mg = 
m-f.i = nif) from 500 GeV to 80 GeV in the MC. For 
niq = nig = 300 GeV MC events, the percentage eo in¬ 
creases by an additional 25%. Sleptons with mass below 
80 GeV have already been excluded [|o| . The increase in 
X 2 production increases the mass limit by approximately 
10 GeV. 

A light stop t\ would also modify squark and gluino 
decays and would therefore affect X 2 production. If m t - 
is lowered from 500 GeV to the lower experimental limit 
of 80 GeV |l(],|ll|], a 15% reduction in X 2 production is 
predicted. This reduction lowers the limit for equal mass 
squarks and gluinos by about 6 GeV. 

Following the above procedure, we obtain a lower 
limit for gluino (squark) mass of 240 GeV when squarks 
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(gluinos) are heavy. Again, these limits vary by approx¬ 
imately 10 GeV if ti and/or sleptons are light. 

In summary, we have searched for an excess of 'y^T 
events with two or more jets in pp collisions at y/s = 
1.8 TeV. Such events are predicted in the minimal su¬ 
persymmetric standard model. We find that the number 
of observed + > 2 jets events agrees well with that 
expected from background processes. Within the frame¬ 
work of the MSSM, with choices of parameters consistent 
with B(x 2 — * 1 2 3 4 5 6 7 > 7X?) = 1 and m^o — m^o > 20 GeV, we 
obtain a 95% C.L. lower mass limit of 310 GeV for equal 
mass squarks and gluinos and of 240 GeV for squarks 
(gluinos) when gluinos (squarks) are heavy. These lim¬ 
its constrain the models discussed in Ref. pj] , but do not 
exclude all of them. 
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